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Abstract 

A significant asymmetry in baryon/antibaryon yields in the central region of high 
energy collisions is observed when the initial state has non-zero baryon charge. This 
asymmetry is connected with the possibility of a baryon charge diffusion in rapidity 
space. Evidently, such a diffusion should decrease the baryon charge in the fragmenta- 
tion region leading to the corresponding decrease of the multiplicity of leading baryons. 
As a result, a new mechanism for Feynman scaling violation in the fragmentation region 
is obtained. We present the quantitative predictions for the Feynman scaling violation 
at LHC energies and even at higher energies that can be important for cosmic ray 
physics. 
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1 Introduction 

The problem of Feynman scaling violation has evident both theoretical and practical 
interest. In particular, this question is very important [H |2] for cosmic ray physics, 
where the difference from the primary radiation to the events registrated on the ground 
or mountain level is determined by the multiple interactions of the so-called leading 
particles (mainly baryons) in the atmosphere. 

Despite the lack of direct measurements of Feynman scaling violation for secondary 
baryon spectra in nucleon-nucleon collisions at energies higher than those of ISR, some 
experimental information from cosmic ray experiments seems to confirm [31 lU [S] the 
presence of significant Feynman scaling violation effects. Now the LHCf Collaboration 
has started the search |6l [7] of Feynman scaling violation effects for the spectra of 
photons (tt^) and neutrons in the fragmentation region at LHC energies. 

In principale, the violation of Feynman scaling in the fragmentation region should 
exist due to the energy conservation, since the spectra of charged particles increase in 
the central region. However, no quantitative predictions can be made without some 
model of the particle production. Actually, even a small decrease of the spectrum of 
some secondaries in a narrow region near xp = 1 may be enough to satisfy the energy 
conservation law. 

The predictions of the Additive Quark Model [S] for Feynman scaling violation in the 
fragmentation region were considered in [9]. In this model, the scaling violation effects 
are connected with the increase of the interaction cross sections, which leads to the 
decrease in the number of quark-spectators which form the spectra of fast secondaries. 
However, to make a description of the energy dependences of the spectra as a function 
of xp additional assumptions and parameters are needed. 

The Quark-Gluon String Model (QGSM) ^U\ has some important advantages in this 
respect, since it is an analytical model and it allows the calculation of the spectra of 
secondaries at different initial energies in the whole xp region. The QGSM is based on 
Dual Topological Unitarization (DTU), Regge phenomenology, and nonperturbative 
notions of QCD. This model is successfully used for the description of multiparticle 
production processes in hadron-hadron [HI [T21 [131 [E], hadron- nucleus [IS [TB], and 
nucleus-nucleus ^7\ collisions. 

In the QGSM high energy interactions are considered as proceeding via the ex- 
change of one or several Pomerons, and all elastic and inelastic processes result from 
cutting through or between Pomerons [18]. Inclusive spectra of hadrons are related 
to the corresponding fragmentation functions of quarks and diquarks, which are con- 
structed using the Reggeon counting rules [19]. The quantitative predictions of the 



QGSM depend on several parameters which were fixed by comparison of the theo- 
retical calculations to the experimental data obtained at fixed target energies. The 
first experimental data obtained at LHC show [20j that the model predictions are in 
reasonable agreement with the data. 

In the frame of QGSM several reasons of Feynman scaling violation in the fragmen- 
tation region exist. The first one is the increase of the average number of exchanged 
Pomerons with the energy, which leads to the corresponding increase of the yields of 
hadron secondaries in the central region and to their decrease in the fragmentation 
region. This effect was considered in [2T1 122]- 

In the case of nuclear (air) targets, the growth of the hN cross section with energy 
leads to the increase of the average number of fast hadron inelastic collisions inside the 
nucleus. Thus, the average number of Pomerons is additionally increased, resulting in 
a stronger Feynman scaling violation [2T| [22] . 

In reference [23] these predictions were taken into account to calculate the penetra- 
tion of fast hadrons into the atmosphere, leading to a better description of the cosmic 
ray experimental data. 

The differences in the yields of baryons and antibaryons produced in the central 
(midrapidity) region of high energy pp interactions [m |20l |2ll |25l |26l [271 [2H] are 
significant. Evidently, the appearance of the positive baryon charge in the central region 
oipp collisions should be compensated by the decrease of the baryon multiplicities in the 
fragmentation region that leads to an additional reason for Feynman scaling violation. 

In the present paper we consider the effects of Feynman scaling violation in pp 
collisions through large distance baryon diffusion in rapidity space. The role of the 
nuclear factor for air nuclei should be of the same magnitude as for that presented in 
Ell [22]. 



2 Baryon/antibaryon asymmetry in the QGSM 

The Quark-Gluon String Model (QGSM) [TOl [HI [12] allows us to make quantitative 
predictions of different features of multiparticle production. In QGSM the inclusive 
spectrum of a secondary hadron h is determined by the convolution of the diquark, 
valence quark, and sea quark distributions u{x,n) in the incident particles with the 
fragmentation functions G^{z) of quarks and di quarks into the secondary hadron h. 
These distributions, as well as the fragmentation functions, are constructed using the 
Reggeon counting rules pS]. The details of the model are presented in (TUl HH [121 Hi]- 
The Pomeron parameters were taken from 



In the string models, baryons are considered as configurations consisting of tfiree 
connected strings (related to three valence quarks) called string junction (SJ) [221 EOl 
31} 132] . Such a baryon structure is supported by lattice calculations [33]. This picture 
leads to some general phenomenological predictions. In the case of inclusive reactions 
the baryon number transfer to large rapidity distances in hadron-nucleon and hadron- 
nucleus reactions can be explained [HI [211 EH EH ES] by S J diffusion. 

The production of a baryon-antibaryon pair in the central region usually occurs 
via SJ-SJ (SJ has upper color indices whereas SJ has lower indices) pair production, 
which then combines with sea quarks and sea antiquarks into a BB pair [31l |37] , as it 
is shown in Fig. la. 
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Figure 1: QGSM diagrams describing secondary baryon B production by diquark d. (a) Central 
production of BB pair. Single B production in the processes of diquark fragmentation: (b) initial SJ 
together with two valence quarks and one sea quark, (c) initial SJ together with one valence quark 
and two sea quarks, and (d) initial SJ together with three sea quarks. Quarks are shown by solid 
curves and SJ by dashed curves. 

In the processes with incident baryons, e.g. in pp collisions, another possibility to 
produce a secondary baryon in the central region exists. This possibility is the diffusion 
in rapidity space of any SJ existing in the initial state and it can lead to significant 
differences in the yields of baryons and antibaryons in the midrapidity region even 
at rather high energies [H] . The most important experimental fact in favour of this 
process is the rather large asymmetry in Q and Q baryon production in high energy 
-K^p interactions [38] . 

The theoretical quantitative description of the baryon number transfer via SJ mech- 



anism was suggested in the 90's and used to predict [2S] the p/p asymmetry at HERA 
energies. 

In order to obtain the net baryon charge we consider, following ref. [H] three 
different possibilities. The first one is the fragmentation of the diquark giving rise to a 
leading baryon (Fig. lb). A second possibility is to produce a leading meson in the first 
break-up of the string and a baryon in a subsequent break-up (Fig. Ic). In these two 
first cases the baryon number transfer is possible only for short distances in rapidity. 
In the third case, shown in Fig. Id, both initial valence quarks recombine with sea 
antiquarks into mesons M while a secondary baryon is formed by the SJ together with 
three sea quarks. 

The fragmentation functions for the secondary baryon B production corresponding 
to the three processes shown in Figs, lb, Ic, and Id, can be written as follows [14j : 

Gl{z) = aM-vf^-z'-', (1) 

G%{z) = aM-vf,-z'.{l-z), (2) 

Giiz) = a^-e-vf,-z'-'^^^-{l-zY, (3) 

where ajy is the normalization parameter, and v^ , f^, f£ are the relative probabilities 
for different baryons production that can be found by simple quark combinatorics 
[lOl HI] . Their numerical values for different secondary baryons were presented in [27] . 

The first two processes shown in Figs, lb and Ic, Eqs. (1) and (2), determine the 
spectra of leading baryons in the fragmentation region. The third contribution shown 
in Fig. Id, Eq. (3), is essential if the value of the intercept of the SJ exchange Regge- 
trajectory, asj, is not too small. In QGSM the weight of this third contribution is 
determined by the coefficient e which fixes the small probability for such a baryon 
number transfer to occur. 

At high energies, the SJ contribution to the inclusive cross section of secondary 
baryon production at large rapidity distance Ay from the incident nucleon can be 
estimated as 

{l/a)da^/dy ^ as ■ e ■ e^^""^-^)-^^ , (4) 

where a^ = ajv ■ w^. The baryon charge transferred to large rapidity distances can be 
determined by integration of Eq. (4), so it is of the order of 

(--) - -B ■ j^^^ , (5) 

and so, only the left part of the initial baryon charge can be used for the production 
of the leading baryons. 



To obtain the QGSM predictions for the spectra of leading baryons we use the 
standard expressions of Reggeon theory and QGSM fTUl \T^ . 

Though currently the value of asj = 0.5 seems more plausible fl2], the value of 
<^SJ = 0.9 can not be excluded [131111]. Thus, in this paper we present the calculation 
obtained with these two values of asj, and also without any SJ contribution [e = 0.). 

3 Spectra of baryons in ep collisions 

The only very high energy reaction where the leading secondary baryons were measured 
is ep collisions at HERA. Here, both ep — )■ epX (i.e. ■yp — )■ pX) [l5l |16] and ep — )■ enX 
(i.e. 7p — )■ nX) [47j were investigated at VT ~ 200 GeV. Due to the vector dominance 
principle, the inelastic 'jp interaction can be considered as a superposition of p^p and 
copu intercations. In the frame of the QGSM these two interactions are equivalent to 
the sum of (l/2)7r+p+ (l/2)7r~p collisions. 

In Fig. 2 the experimental data for 'jp — > pX ^5\ H6] and for 7p — )■ nX [17] are 
compared with the QGSM calculations. Here the experimental data are presented 
as depending on xl which is very close to xp if Pt is small. The comparison of the 
ep — )■ epX data with the QGSM calculations is shown in the upper panel of Fig. 2, 
while the data of the ep — )■ enX reaction at low Q^ and high Q^ are presented in the 
middle and lower panels of Fig. 2, respectively. 

The solid curves correspond to the calculations with asj = 0.9, dashed curves 
to those with asj = 0.5, and dotted curves are calculated without SJ contributions 
{e = 0.). Sice the energy dependences of all curves are weak, they all are calculated at 
W = 200 GeV. The difference between the curves is small. All experimental points are 
obtained at rather small values of the transverse momenta of the secondary baryons, 
so they lie below the theoretical curves that are calculated for the spectra integrated 
over Pt- 

Usually, the spectra of neutrons produced in DIS at large xl {xl > 0.7) and small 
Pt are described in the framework of the one-pion-exchange approach [IHl SS]. This 
approach leads to a more detailed description of the data, for example it allows one to 
calculate the x^-spectra in different pt regions. However, our QGSM description based 
on quark and diquark fragmentation picture gives reasonable results for xl < 0.8. 



^As we do not consider the yields of strange secondaries, the numerically small contribution of ss 
pairs (0-meson) is not important here. 
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Figure 2: The QGSM predictions for the spectra of secondary protons (upper panel) and neutrons 
(middle and lower panels) produced in jp collisions, together with the experimental data of references 
[ISlll^ and [47]. Solid curves correspond to the value asj — 0.9, dashed curves to the value asj — 0.5, 
and dotted curves are calculated without SJ contribution (e = 0.). 



4 Predictions for the spectra of baryons in pp colli- 



sions 



The QGSM predictions for the inclusive spectra of secondary protons, neutrons, and 
A at energies v^ = 200 GeV, 900 GeV, 7 TeV, and 100 TeV are presented in Figs. 3, 
4, and 5. In all the calculations we have accounted for the exact conservation of the 
baryon charge. 

In all Figures 3, 4, and 5 one can see the peaks at very small xp, that are connected 
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Figure 3: The QGSM predictions for the spectra of secondary protons at energies ^/s = 200 GeV (up 
left panel), ^/s = 900 GeV (up right panel), ^/s = 7 TeV (down left panel), and ^i = 100 TeV (down 
right panel). Solid curves correspond to the value asj ~ 0.9, dashed curves to the value asj — 0.5, 
and dotted curves are calculated without SJ contribution (e = 0.). 



to the contribution of the pair BB production via the mechanism shown in Fig. la. 
The increase of the spectra with xp is connected to the contributions of the processes 
shown in Figs, lb and Ic. The difference between both the sohd and dashed curves 
to the dotted curves show the effect of baryon number transfer to small Xp region due 
to the diagram Fig. Id. An important result is that the Feynman scaling violation is 
more sensitive to the fact of the inclusion of the baryon number diffusion than to the 
exact value of asj- 

One can see that in the case of secondary protons, shown in Fig. 3, the differences 
in the fragmentation region between the two calculations with SJ contribution, with 
values asj = 0.9 and asj = 0.5, are very small at energies y/s < 1 TeV, while they 
increase up to 5-7% at higher energies. In the case of the calculations with and without 
SJ contribution the difference in the spectra in the region xp = 0.5 — 0.8 is of about 
7-10% for yi = 200 GeV and 900 GeV, and it becomes of about 15-20% at v^ = 7 
TeV and 100 TeV. So, by accounting for the possibility of baryon charge diffusion in 
rapidity space one gets one additional contribution to the Fyenman scaling violation. 
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Figure 4: The QGSM predictions for the spectra of secondary neutrons at energies y/s = 200 GeV 
(up left panel), -^/s = 900 GcV (up right panel), y/s ~ 7 TcV (down left panel), and y/s = 100 
TcV (down right panel). Solid curves correspond to the value asj — 0.9, dashed curves to the value 
as J — 0-5, and dotted curves are calculated without SJ contribution (e = 0.). 



Our predictions for the spectra of secondary neutrons are presented in Fig. 4. Their 
behaviour with respect to the accounting of the SJ contribution and to the effect of 
the Feynman scahng violation are similar to those of secondary protons. 

In the case of secondary A-hyperon production, shown in Fig. 5, the spectra de- 
crease rather fast at large xp due to the faster decrease of uu and ud fragmentation 
functions into A in comparison with their fragmentation into secondary nucleon, and 
also to the relatively smaller contribution of diffraction dissociation. The absolute val- 
ues of A spectra are smaller than those of the protons and neutrons spectra due to the 
strangeness suppression factor. 

The Feynman scaling violation effects are shown in more detail in Figs. 6 and 7, for 
secondary protons and neutrons, respectively. Here, we present the energy dependences 
of the spectra at four values of xp, namely at xp = 0.05, 0.2, 0.5, and 0.7. Separately, 
we present for the first two values of xp (i.e. in the top panels of Figs. 6 and 7) the 
energy dependences of net baryon production, i.e. the differences p — p and n — n. 
In the case of xp = 0.5 and 0.7 the differences between the baryon spectra and the 
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Figure 5: The QGSM predictions for the spectra of secondary A-hyperons at energies ^/s = 200 
GeV (up left panel), y/s — 900 GeV (up right panel), y/s = 7 TeV (down left panel), and ^/s — 100 
TeV (down right panel). Solid curves correspond to the value asj — 0.9, dashed curves to the value 
as J — 0.5, and dotted curves are calculated without SJ contribution (e — 0.). 

net-baryon spectra are negligible. 

One can see that at xp = 0.05 and 0.2 the spectra of secondary protons and neu- 
trons increase with energy (both the total spectra as well as the net baryon spectra). 
At larger xp these spectra decrease with the energy. These energy dependences are 
connected with the growth of the average number of exchanging Pomerons. The differ- 
ences between solid and dashed curves to dotted curves show the effect of SJ diffusion. 



5 The ratios of the inclusive spectra of baryons at 
different energies 

The possible Feynman scaling violation at superhigh energies was discussed in [6], 
based on Monte Carlo calculations. These effects are, as a rule, numerically not large, 
and so it is more suitable to consider the ratios of the spectra at different energies. 
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Figure 6: The QGSM predictions for the spectra of secondary protons as the functions of energy 
at fixed values of xp- Thin curves show the total proton spectra and bold curves at top panels 
the spectra of net protons, i.e. the values of the p — p differences. Solid curves correspond to the 
value asj — 0.9, dashed curves to the value asj — 0.5, and dotted curves are calculated without SJ 
contribution (e = 0.). 



These ratios calculated in the QGSM at y/s = 7 TeV, 14 TeV, and 1000 TeV to the 
values at y/s = 900 GeV, are presented both for secondary protons and neutrons in 
Fig. 8. 

In Fig. 9 we present the ratios of the spectra of secondary protons and neutrons 
as the functions of rapidity, calculated with asj = 0.9 and without SJ contribution, 
at ^/s = 200 GeV (solid curves) and at ^/s = 100 TeV (dashed curves). The similar 
ratios calculated with asj = 0.5 and without SJ contribution are shown by dash-dotted 
curves at ^/s = 200 GeV and by dotted curves at ^/s = 100 TeV. 



6 Conclusion 

We present the QGSM predictions for Feynman scaling violation in the spectra of 
leading baryons due to baryon charge diffusion at large distances in the rapidity space. 
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Figure 7: The QGSM predictions for the spectra of secondary neutrons as the functions of energy 
at fixed values of xp- Thin curves show the total neutron spectra and bold curves at top panels 
the spectra of net neutrons, i.e. the values of n — n differences. Solid curves correspond to the 
valuea^j — 0.9, dashed curves to the value asj — 0.5, and dotted curves are calculated without SJ 
contribution (e — 0.). 

The existance of such a diffusion has been attested observed in many papers, even at 
the LHC energies, and so the decrease of the spectra in the fragmentation region is the 
direct consequence of baryon charge conservation. However, the numerical values of 
the scaling violation at different Xp are model-dependent. On the other point of view, 
the calculations of scaling violation for leading baryons should be accompanied by the 
calculation of the baryon/antibaryon asymmetry in the central region. 

The first experimental data obtained at LHC are in general agreement with the 
QGSM calculations performed with the same values of parameters which were deter- 
mined at lower energies (mainly for the description of the the fixed target experiments). 
However, the numerical value of agj is not well-known. ALICE Collaboration data on 
the ratios of pp production are in agreement with a value asj = 0.5, while LHCb 
Collaboration data on the ratios of AA production do not allow to fix the value of asj-, 
and there are two experimental results for BB production asymmetry. 

We neglect by the possibility of interactions between Pomerons the (so-called en- 
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Figure 8: The QGSM predictions for the ratios of the spectra of secondary protons (left panels) and 
neutrons (right panels) at three different energies to those at y^ = 900 GeV. 

hancement diagrams), since our estimations [50] show that the inclusive density of 
secondaries produced in pp colhsions at LHC energies is not large enough for these 
diagrams to be significant. 

Our calculations are in reasonable agreement with the results of ref. 
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Figure 9: The QGSM predictions for the ratios of the spectra of secondary protons (left panels) 
and neutrons (right panels) calculated with asj = 0.9 (solid and dashed curves) and with asj = 0.9 
(dash-dotted and dotted curves) and with energies -^i = 200 GeV (solid and dash-dotted curves) and 
at ^/s = 100 TeV (dashed and dotted curves). 
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